12 January 2013 



Elliptic Flow: A Brief Review 

Raimond Snellings 

Utrecht University, P.O. Box 80000, 3508 TA Utrecht, The Netherlands 

Abstract. One of the fundamental questions in the field of subatomic physics is what 
happens to matter at extreme densities and temperatures as may have existed in the 
first microseconds after the Big Bang and exists, perhaps, in the core of dense neutron 
stars. The aim of heavy-ion physics is to collide nuclei at very high energies and thereby 
create such a state of matter in the laboratory. The experimental program started 
in the 1990's with collisions made available at the Brookhaven Alternating Gradient 
Synchrotron (AGS), the CERN Super Proton Synchrotron (SPS) and continued at 
the Brookhaven Relativistic Heavy-Ion Collider (RHIC) with maximum center of mass 
energies of ^/s NN = 4.75, 17.2 and 200 GeV respectively. Collisions of heavy-ions at the 
unprecedented energy of 2.76 TeV have recently become available at the LHC collider at 
CERN. In this review I will give a brief introduction to the physics of ultra-relativistic 
heavy-ion collisions and discuss the current status of elliptic flow measurements. 
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1. Heavy-Ion Physics 

To our current understanding, the universe went through a series of phase transitions 
which mark the most important epochs of the expanding universe after the Big Bang. 
At 10~ n s and at a temperature of T ~ 100 GeV (~ 10 15 K) the electroweak phase 
transition took place where most of the known elementary particles acquired their Higgs 
masses [HOE]. At 10^ 5 s and at a temperature of ~ 200 MeV (~ 10 12 K), the strong 
phase transition took place where the quarks and gluons became confined into hadrons 
and where the approximate chiral symmetry was spontaneously broken [1]. 

Quantum Chromodynamics (QCD) is the underlying theory of the strong force. 
Although its fundamental degrees of freedom (quarks and gluons) cannot be observed as 
free particles the QCD Lagrangian is well established. One of the key features of QCD 
is the self coupling of the gauge bosons (gluons) which cause the coupling constant to 
increase with decreasing momentum transfer. This running of the coupling constant 
gives rise to asymptotic freedom [SJ E] and confinement at large and small momentum 
transfers, respectively. At small momentum transfer nonperturbative corrections, which 
are notoriously hard to calculate, become important. For this reason two important 
nonperturbative properties of QCD, confinement and chiral symmetry breaking, are 
still poorly understood from first principles. 

One of the fundamental questions in QCD phenomenology is what the properties of 
matter are at the extreme densities and temperatures where the quarks and gluons are 
in a deconfined state, the so-called Quark Gluon Plasma (QGP) |7J. Basic arguments [7] 
allow us to estimate the energy density e ~ 1 GeV/fm 3 and temperature T ~ 200 MeV 
at which the strong phase transition takes place. These values imply that the transition 
occurs in a regime where the coupling constant is large so that we can not rely anymore 
on perturbative QCD. Better understanding of the non-perturbative domain comes from 
lattice QCD, where the field equations are solved numerically on a discrete space-time 
grid. Lattice QCD provides quantitative information on the QCD phase transition and 
the Equation of State (EoS) of the deconfined state. At extreme temperatures (large 
momenta) we expect that the quarks and gluons are weakly interacting and that the 
QGP would behave as an ideal gas. For an ideal massless gas the EoS is given by: 

1 7T 2 

P=-e, e = g—T\ (1) 
3 ' y 30 ' K J 

where P is the pressure, e the energy density, T the temperature and g is the effective 

number of degrees of freedom. Each bosonic degree of freedom contributes 1 unit to g, 

whereas each fermionic degree of freedom contributes |. The value of g = 47.5 for a 

three flavor QGP which is an order of magnitude larger than that of a pion gas where 

g~ 3. 

Figure [T] shows the temperature dependence of the energy density as calculated from 
lattice QCD [8J. It is seen that the energy density changes rapidly around T ~ 190 MeV, 
which is due to the rapid increase in the effective degrees of freedom (lattice calculations 
show that the transition is a crossover). Also shown in Fig. [l] is the pressure which 
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Figure 1. Energy density e/T 4 (full curve) and pressure 3P/T 4 (dashed curve) as a 
function of temperature T from lattice calculations [5] . The arrow indicates the Stefan 
Boltzmann limit of the energy density. 

changes slowly compared to the rapid increase of the energy density around T = 190 
MeV. It follows that the speed of sound, c s = ^dP/de, is reduced during the strong 
phase transition. At large temperature the energy density reaches a significant fraction 
(~ 0.9) of the ideal massless gas limit (Stefan-Boltzmann limit). 

Relativistic heavy-ion collisions are a unique tool to create and study hot QCD 
matter and its phase transition under controlled conditions [U El QUI HH [12], [13]. As 
in the early universe, the hot and dense system created in a heavy-ion collision will 
expand and cool down. During this evolution the system probes a range of energy 
densities and temperatures, and possibly different phases. Provided that the quarks 
and gluons undergo multiple interactions the system will thermalize and form the QGP 
which subsequently undergoes a collective expansion and eventually becomes so dilute 
that it hadronizes. This collective expansion is called flow. 

Flow is an observable that provides experimental information on the equation of 
state and the transport properties of the created QGP. The azimuthal anisotropy in 
particle production is the clearest experimental signature of collective flow in heavy-ion 
collisions [HJ [151 US El HE] • This so-called anisotropic flow is caused by the initial 
asymmetries in the geometry of the system produced in a non-central collision. The 
second Fourier coefficient of the azimuthal asymmetry is called elliptic flow. In this 
report I will describe the relation between elliptic flow and the geometry of the collision 
(section [2] and [3]) , the sensitivity of elliptic flow to the EoS and transport properties 
(section [3]) and the techniques used to measure elliptic flow from the data (section [4]). 
In sectional will review the elliptic flow measurements at the LHC and at lower energies, 
together with the current theoretical understanding of these results. 
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Figure 2. Left: The two heavy-ions before collision with impact parameter b. Right: 
The spectators continue unaffected, while in the participant zone particle production 
takes place. 



Heavy-ions are extended objects and the system created in a head-on collision 
is different from that in a peripheral collision. To study the properties of the 
created system, collisions are therefore categorized by their centrality. Theoretically 
the centrality is defined by the impact parameter b (see Fig. [2]) which, however, 
cannot be directly observed. Experimentally, the collision centrality can be inferred 
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Figure 3. a) Charged particle distribution from Pb-Pb collisions at y/s NN = 2.76 TeV 
measured with ALICE, showing a classification in centrality percentiles (from [20] ) . 
b) Number of participating nucleons -ZV par t and binary collisions A^i n versus impact 



parameter for Pb-Pb and Au-Au collisions at -Vs^ 



2.76 and 0.2 TeV, respectively. 



from the measured particle multiplicities, given the assumption that the multiplicity 
is a monotonic function of b. The centrality is then characterized by the fraction, 
7r6 2 /7r(2i?yi) 2 , of the geometrical cross-section with Ra the nuclear radius (see Fig. pm). 
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Instead of by impact parameter, the centrality is also often characterized by the 
number of participating nucleons (nucleons that undergo at least one inelastic collision) 
or by the number of equivalent binary collisions. Phenomenologically it is found that the 
total particle production scales with the number of participating nucleons whereas hard 
processes scale with the number of binary collisions. These measures can be related 
to the impact parameter b using a realistic description of the nuclear geometry in a 
Glauber calculation [Tj5], as is shown in Fig. [3]d. This Figure also shows that Pb-Pb 
collisions at y/s NN =2.76 TeV and Au-Au at y/s NN = 0.2 TeV have a similar distribution 
of participating nucleons. The number of binary collisions increases from Au-Au to Pb- 
Pb by about 50% because the nucleon-nucleon inelastic cross section increases by about 
that amount at the respective center of mass energies of 0.2 and 2.76 TeV. 

3. Anisotropic Flow 

Flow signals the presence of multiple interactions between the constituents of the 
medium created in the collision. More interactions usually leads to a larger magnitude 
of the flow and brings the system closer to thermalization. The magnitude of the 
flow is therefore a detailed probe of the level of thermalization. The theoretical tools 




Figure 4. Almond shaped interaction volume after a non-central collision of two 
nuclei. The spatial anisotropy with respect to the x-z plane (reaction plane) translates 
into a momentum anisotropy of the produced particles (anisotropic flow). 

to describe flow are hydrodynamics or microscopic transport (cascade) models. In the 
transport models flow depends on the opacity of the medium, be it partonic or hadronic. 
Hydrodynamics becomes applicable when the mean free path of the particles is much 
smaller than the system size, and allows for a description of the system in terms of 
macroscopic quantities. This gives a handle on the equation of state of the flowing 
matter and, in particular, on the value of the sound velocity c s . 

Experimentally, the most direct evidence of flow comes from the observation of 
anisotropic flow which is the anisotropy in particle momentum distributions correlated 
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Figure 5. The created initial transverse energy density profile and its time 
dependence in coordinate space for a non-central heavy-ion collision |21) . The z-axis 
is along the colliding beams, the x-axis is defined by the impact parameter. 



with the reaction plane. The reaction plane is defined by the impact parameter and the 
beam direction z (see Fig. [4]) . A convenient way of characterizing the various patterns 
of anisotropic flow is to use a Fourier expansion of the invariant triple differential 
distributions: 

d 3 N 1 d 2 iV / ~ \ 

where E is the energy of the particle, p the momentum, p t the transverse momentum, (p 
the azimuthal angle, y the rapidity, and \1/rp the reaction plane angle. The sine terms 
in such an expansion vanish because of the reflection symmetry with respect to the 
reaction plane. The Fourier coefficients are p t and y dependent and are given by 

v n (Pt,y) = (cos[n((p - *rp)]), (3) 

where the angular brackets denote an average over the particles, summed over all events, 
in the (pt,y) bin under study. In this Fourier decomposition, the coefficients V\ and v 2 
are known as directed and elliptic flow, respectively. 

The evolution of the almond shaped interaction volume is shown in Fig. [5j The 
contours indicate the energy density profile and the plots from left to right show how 
the system evolves from an almond shaped transverse overlap region into an almost 
symmetric system. During this expansion, governed by the velocity of sound, the created 
hot and dense system cools down. 

Figure |6^i shows the velocity of sound versus temperature for three different 
equations of state [22J. The dash-dotted line is the hadron resonance gas EoS, the 
red full line is a parameterization of the EoS which matches recent lattice calculations 
and the blue dashed line is an EoS which incorporates a first order phase transition. 
The arrows indicate the corresponding transition temperatures for the lattice inspired 
EoS and the EoS with a first order phase transition. The temperature dependence of 
the sound velocity clearly differs significantly between the different equations of state. 
Because the expansion of the system and the buildup of collective motion depend on the 
velocity of sound, it is expected that this difference will have a clear signature in the flow. 
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Figure 6. a) The velocity of sound squared versus temperature for three equations 
of state [22) . b) The anisotropy in momentum space for two equations of state used in 
hydrodynamic calculations |22j . 



The buildup of the flow for two different EoS is shown in Fig. [6}d. Due to the stronger 
expansion in the reaction plane the initial almond shape anisotropy in coordinate space 
vanishes, as was shown in Fig. |5j while the momentum space distribution changes in 
the opposite direction from being approximately azimuthally symmetric to having a 
preferred direction in the reaction plane. The asymmetry in momentum space can be 
quantified by: 

£ P = Jf +T y l 4 J 

\J-xx T J-yy/ 

where T xx and T yy are the diagonal transverse components of the energy momentum 
tensor and the brackets denote an averaging over the transverse plane. Figure [6]d shows 
that e p versus time starts at zero after which the anisotropy quickly develops and is 
indeed dependent on the EoS. 

Although e p is not a direct observable, the observed EoS dependence of e p versus 
time is reflected in the experimental observable i n particular when plotted as function 
of transverse momentum and particle mass. Figure [7^i shows ^-differential elliptic flow 
for pions and protons after the transverse momentum spectra have been constrained. 
A clear mass dependence of v 2 at low transverse momentum is observed for both 
equations of state. The figure also clearly shows that the pion v 2 does not change 
much between the lattice EoS and EoS Q. On the other hand, the v 2 of protons does 
change significantly because the heavier particles are more sensitive to the change in 
collective motion. Therefore measurements of V2(pt) for various particle species provide 
an excellent constraint on the EoS in ideal hydrodynamics. 

More recently, it was realized that small deviations from ideal hydrodynamics, in 
particular viscous corrections, already modify significantly the buildup of the elliptic 
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Figure 7. a) The EoS dependence of V2{pt) for pions and protons [22]. The 
full lines are for the lattice inspired EoS and the dashed lines for an EoS which 
incorporates a first order phase transition, b) The dependence on rj/s of V2(jh) fo r 
charged particles [24] . 



flow [23]. The shear viscosity determines how good a fluid i^JJ however, for relativistic 
fluids the more useful quantity is the shear viscosity over entropy ratio rj/s. Known good 
fluids in nature have an rj/s of order h/hs- In a strongly coupled Af = 4 supersymmetric 
Yang Mills theory with a large number of colors ('t Hooft limit), r]/s can be calculated 
using a gauge gravity duality [25J: 

r] h 

s Ank-Q 

Kovtun, Son and Starinets conjectured, using the AdS/CFT correspondence, that this 
implies that all fluids have rj/s > Ti/Ank-Q (the KSS bound.). We therefore call a fluid 
with rj/s — 1/Att (in natural units) a perfect fluid. The KSS bound raises the interesting 
question on how fundamental this value is in nature and if the QGP behaves like an 
almost perfect fluid. It is argued that the transition from hadrons to quarks and gluons 
occurs in the vicinity of the minimum in rj/s, just as is the case for the phase transitions 
in helium, nitrogen, and water. An experimental measurement of the minimal value of 
rj/s would thus pinpoint the location of the transition [2S1 12"T] . 

Experimentally we might get an answer to the magnitude of rj/s by measuring vi 
as is shown in Fig. The full line is close to ideal hydrodynamics (rj/s ~ 0) while 
the three other lines correspond to rj/s values of up to three times the KSS bound. 
Different magnitudes of rj/s clearly lead to a dramatically different magnitude of t>2 and 
change its p t dependence. However the magnitude and p t dependence of v 2 not only 
depend on rj/s but also on the EoS as we have seen. 

The magnitude of v% does not only depend on the medium properties of interest, 
but is also proportional to the initial spatial anisotropy of the collision region. This 

\ a good fluid has a small viscosity and does not convert much kinetic energy of the flow into heat. 
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Figure 8. a) The eccentricity e calculated in a color glass condensate (CGC) model 
and using a Glauber model [28] , b) The v-i obtained using the CGC or Glauber initial 
eccentricity [28] . 



spatial anisotropy can be characterized by the eccentricity, which is defined by 

(y 2 - x 2 ) 



(5) 



(y 2 + x 2 ) ' 

where x and y are the positions of the participating nucleons in the transverse plane 
and the brackets denote an average which traditionally was taken over the number of 
participants. Recent calculations have shown that the eccentricity obtained in different 
descriptions, in particular comparing a Glauber with a Color Glass Condensate (CGC) 
description, shows that e varies by almost 25% at a given impact parameter [28J, see 
Fig. |8^l. The elliptic flow, obtained when using these different initial eccentricities is 
shown in Fig. ^jp. As expected, the different magnitude of the eccentricity propagates to 
the magnitude of the elliptic flow. Because currently we cannot measure the eccentricity 
independently this leads to a large uncertainty in experimental determination of 77/s. 

To summarize, we have seen that the elliptic flow depends on fundamental 
properties of the created matter, in particular the sound velocity and the shear viscosity, 
but also on the initial spatial eccentricity. Detailed measurements of elliptic flow as 
function of transverse momentum, particle mass and collision centrality provide an 
experimental handle on these properties. In the next section, before we discuss the 
measurements, we first explain how we estimate the anisotropic flow experimentally. 



4. Elliptic Flow: Analysis Methods 

Because the reaction plane angle is not a direct observable the elliptic flow (Eq. [3]) can 
not be measured directly so that it is usually estimated using azimuthal correlations 
between the observed particles. Two-particle azimuthal correlations, for example, can 
be written as: 
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= (( e » 2 (^-*Rp))( e -* 2 (v2-* R p)) + ( j 2 ) ) 

= <«| + fc>, (6) 

where the double brackets denote an average over all particles within an event, followed 
by averaging over all events. In Eq. [6] we have factorized the azimuthal correlation 
between the particles in a common correlation with the reaction plane (elliptic flow v^) 
and a correlation independent of the reaction plane (non-flow 62). Here we have assumed 
that the correlation between v 2 and 62 is negligible. If 62 is small, Eq. [6] can be used 
to measure but in general the non-flow contribution is not negligible. In Fig. M 






Figure 9. Examples of particle distributions in the transverse plane, where for a) 
v 2 > 0, v 2 {2} > 0, b) v 2 = 0, v 2 {2} = 0, and c) v 2 = 0, v 2 {2} > 0. 

we illustrate two-particle nonflow contributions as follows: In Fig. [9^i an anisotropic 
distribution is shown for which both v% = (cos 20) and the two-particle correlation 
v 2 {2} = \J (cos 2(0x — 2 )) are positive. Figure 9Jd shows a symmetric distribution for 
which v 2 = and also t>2{2} = 0. Figure]^ shows two symmetric distributions rotated 
with respect to each other which give V2 = while v 2 {2} is nonzero. This illustrates 
how non-flow contributions from sources like resonance decays or jets can contribute to 
t>2 measured from two particle correlations. 

The collective nature of elliptic flow can be exploited to suppress non-flow 
contributions [291 EQ]. This is done using so called cumulants, which are genuine multi- 
particle correlations. For instance, the two particle cumulant C2{2} and the four particle 
cumulants C2{4} are defined as: 



c 2 {2} = 
c 2 {4} = 



,i2(i^l+«p2— ¥>3— <pa) 



(7) 



V2+64 + 4«2<J 2 + 2S 2 



v\ + 5 2 



4 + h 



(8) 

From the combinatorics it is easy to show that 82 oc 1/M C and £4 oc where M c is 

the number of independent particle clusters. Therefore, t>2{2} is only a good estimate 
if V2 3> \ j \[W C while t>2{4} is already a good estimate of v 2 if v 2 3> 1/M C 3//4 ; for c 2 {oo} 
this argument leads to t> 2 3> 1/M C . This shows that for a typical Pb-Pb collision at 
the LHC with M c = 500 the possible non-flow contribution can be reduced by more 
than an order of magnitude using higher order cumulants. One of the problems in using 
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multi-particle correlations is the computing power needed to go over all possible particle 
multiplets. To avoid this problem, multi-particle correlations in heavy-ion collision are 
calculated from generating functions with numerical interpolations [29] or, as was shown 
more recently, from an exact solution [31 J. 

The last equality in Eq. [8] follows from the assumption that v 2 and 5 2 are 
uncorrelated and also that (6 2 ) = (i^) 2 and (v 2 ) = (v 2 ) 2 . In other words, we have 
neglected the event- by-event fluctuations in v 2 and 5 2 . The effect of the fluctuations on 
v 2 estimates can be obtained from 

(v 2 2 ) = (v 2 ) 2 + a 2 , 

(vj) = (v 2 ) A + 6a 2 (v 2 ) 2 , 

(v e 2 ) = (v 2 ) 6 + 15a 2 (v 2 ) A . (9) 
Neglecting the non-flow terms we have the following expressions for the cumulants: 

M2} = 



«a{4} = p{v 2 2 ? - (n 



v 2 {6} = ^\ ((vl) - 9(vl)(vl) + 12(^)3). (10) 

Here we have introduced the notation v 2 {n} as the flow estimate from the cumulant 
ing that cr <^ (i 

v 2 {2} = (v 2 ) + 



c 2 {n}. Assuming that cr <C (v) we obtain from Eqs. [9] and 10, up to order a 2 : 

1 a 2 



M4} = (v 2 ) ~ 



2W 
1 a 2 



2(v 2 ) 



MG} = (V 2 )-I£y (ID 



From Eqs. [7] and 11 it is clear that the difference between v 2 {2} and t>2{4} is sensitive 
to non-flow and fluctuations. 

Flow fluctuations have become an important part of elliptic flow studies 



EH EHl [361 [371 EHl [391 SOI SH S2]. It is believed that such fluctuations originate mostly 
from fluctuations in the initial collision geometry. This is illustrated in Fig. 10 which 
shows participants that are randomly distributed in the overlap region. This collection 
of participants defines a participant plane ^pp [33J which fluctuates, for each event, 
around the reaction plane \I/rp. These fluctuations can be estimated from calculations 
in, for instance, a Glauber model. 



Figure Hi shows the eccentricities (Eq. [5J) calculated in a Glauber model. Here 
e{RP} denotes the eccentricity in the reaction plane, e is the participant eccentricity 
and e{2} and e{4} are the participant eccentricities calculated using the cumulants, 
analogous to the definitions in Eq. 10 [32J. In Fig. 11 1 the eccentricities are calculated 
using as a weight the participating nucleons (open and solid markers) or as a weight 
binary collisions (dashed lines). The figure clearly shows that in both cases e is in 
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Figure 10. Transverse view of a heavy-ion collision with the reaction plane V^rp 
oriented along the x-axis. Indicated are the participants in the overlap region that 
randomly define a particpant plane ^pp for each collision. 




centrality percentile centrality percentile 

Figure 11. a) The eccentricities from a Glauber calculation for participating nuclcons 
(the solid and open markers) and binary collisions (the dashed lines), b) Various i>2 
estimates compared to the reaction plane value, «2{RP}. 



between e{2} and e{4} as is expected from Eq. 11 The figure also shows that e{4} 
is close to e{RP} for the — 40% centrality range [IUJ HI]. In Fig. 11 d we show a 
transport model calculation of t>2 in the AMPT model [IB]. In this model the true 



reaction plane is known so that we can compare the different v<i estimates with the 
value in the reaction plane. The AMPT model uses a Glauber model for the initial 
conditions and we can therefore compare these estimates with Fig. 11 1 (the dashed lines). 
The agreement between t>2{4} and ^{RP} holds for most of the centrality range, while 
for the eccentricities in the Glauber model a large difference is observed for the more 
peripheral collisions 
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5. Elliptic Flow Measurements 

5.1. The Perfect Liquid 

The large elliptic flow observed at RHIC provides compelling evidence for strongly 
interacting matter which appears to behave like an almost perfect liquid (HI US]. To 
quantify the agreement with an almost perfect fluid the significant viscous corrections 
need to be calculated. Based on different model assumptions the ratio rj/s has 
been estimated at y/s NN = 200 GeV and is found to be below five times the KSS 
bound [TH SZl EH [S]. 




(l/S overia /dN ch /dY)[fm 2 ] P t [GeV/c] 



Figure 12. a) The centrality dependence of t>2{2} compared to viscous hydrodynamic 
model calculations 24 . b) The transverse momentum dependence of v-i compared to 
the same viscous hydrodynamic calculations 24J. 



In Fig. 12 we show the centrality and transverse momentum dependence of t>2 
compared to viscous hydrodynamic calculations |23] with different values of rj/s. Using 
an eccentricity from a CGC inspired calculation, it is seen that both the centrality 
and transverse momentum dependence is well described with an rj/s of two times the 
KSS bound. These calculations are performed under the assumption that the value 
of rj/s is constant during the entire evolution. The value used in these calculations 
should be considered as an effective average of rj/s, because we know from other fluids 
that rj/s depends on temperature. In addition, we also know that part of the elliptic 
flow originates from the hadronic phase. Therefore, knowledge of the temperature 
dependence and knowledge of the relative contributions from the partonic and hadronic 
phase is required to quantify rj/s of the partonic fluid. 

Not only the v 2 of charged particles but also that of identified particles at RHIC 
is described in the framework of viscous hydrodynamics at low-p t . Figure 13 shows the 
measured pion and proton elliptic flow measured by STAR compared to VISHNU [50] 
model calculations. The VISHNU model is a hybrid model which uses viscous 
hydrodynamics for the initial stage followed by a hadron cascade afterburner. In the 
initial viscous hydrodynamic stage rj/s is temperature independent. The rj/s magnitude 
required to describe the pion and proton elliptic flow data are found to be one or two 
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Figure 13. The t>2(j>t) for pions and protons measured by STAR compared to model 
calculations with different eccentricities and rj/s (50] . 



times the KSS bound for a Glauber or CGC eccentricity, respectively. This is in 
agreement with the magnitude of rj/s required to describe charged particle vi- 

While the description of t>2 measurements is encouraging, it is important to 
realize that there are still large uncertainties in: (i) the initial eccentricity, (ii) the 
relative contributions from the hadronic and partonic phase and (iii) the temperature 
dependence of rj/s. Elliptic flow measurements at the LHC, with a higher center of 
mass energy, will constrain these uncertainties and will eventually provide a decisive 
test which of the currently successful model descriptions is the more appropriate. 

5.2. Energy Dependence 

Lead-Lead collisions at the LHC are expected to produce a system which is hotter and 
has a longer lived partonic phase than the system created in Au-Au collisions at RHIC 
energies. As a consequence, the hadronic contribution to the elliptic flow decreases which 
reduces the uncertainty on the determination of rj/s in the partonic fluid. Because rj/s 
is expected to depend on temperature in both the partonic and hadronic system it was 
not clear if the elliptic flow would increase or decrease in going from RHIC to LHC 
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energies. Hydrodynamic models jHU E21 E3] and hybrid models [SJ [55] that successfully 



_£i 0.08 



0.06 - 

0.04 

0.02 '- 
~ 
-0.02 '- 
-0.04 j- 
-0.06 
-0.08 '- 



V 



* ALICE 
TV STAR 

* PHOBOS 
□ PHEN1X 
■ NA49 

O CERES 

+ E877 

* EOS 
A E895 
T FOPI 



10 



10* 



10 J 10 4 



0.12 

0.1 h 
0.08 
0.06 
0.04 
0.02 




~ 10-30% in v 2 . 




,,,,,,,,,,,,,,,,,,,,, 


,,,,,,, 




• • 


• * j 


































• v 2 {2> 






■ v 2 {4} 






Luzum iys=0.08 (Glauber) 






Luzum rys=0.16 (fKLN) 






— Heinz iys=0.20 (fKLN) 




V. . . i . . . . i 


. ... i .... i .... i .... i . 





10 20 30 40 50 60 70 80 
centrality percentile 



Figure 14. a) Integrated elliptic flow at 2.76 TeV in the 20-30% centrality class 
compared with results from lower energies taken at similar centralities (From [10]). b) 
Elliptic flow as a function of event centrality, for the 2- (full circles ) and 4-particle (full 
squares) cumulant methods compared to viscous hydrodynamic calculations (dashed 
lines) [201ISB1IS7]. 



Figure [14p , shows the measured integrated elliptic flow at the LHC in one centrality 
bin, compared to results from lower energies. It shows that there is a continuous increase 
in the elliptic flow from RHIC to LHC energies. In comparison to the elliptic flow 
measurements in Au-Au collisions at y/s NN = 200 GeV v 2 increases by about 30% at 



^ = 2.76 TeV. 



Figure [T4b shows the v 2 for different centralities measured by ALICE with the two- 
and four-particle cumulant method. The difference between the two- and four-particle 
flow estimates for the more central collisions (< 40%) is expected to be dominated by 



event-by-event flow fluctuations (see Eq. 11). For the more peripheral collisions the 
two-particle cumulant is likely biased by non-flow. We already mentioned that t>2{4} 
yields estimates of the elliptic flow in the reaction plane which can thus be compared to 



model predictions of f2{RP}. The curves in Fig. 14 3 show t> 2 {RP} from hydrodynamic 
model calculations for 



s NN = 2.76 TeV, with initial eccentricities and magnitudes of 
rj/s which described the RHIC data. It is seen that in hydrodynamic calculations the 
observed increase in v 2 from RHIC to LHC energies is within expectations. Detailed 
comparisons, however, have to wait till measurements of identified particle spectra and 
identified particle elliptic flow become available. It will then be important to see if one 
still obtains a quantitative description of the data in viscous hydrodynamics and what 
the required magnitude of r]/s then will be. 

In addition to comparisons with detailed dynamic model calculations we might 
also learn something from what happens to the several simple scaling properties 
observed at lower energies. For instance, it was shown that the integrated elliptic flow 
depends linearly on the pseudorapidity rj, measured with respect to the beam rapidity 

1/beam 



as is shown in Fig. 15 1. Based on this scaling behaviour a phenomenological 



Elliptic Flow: A Brief Review 



16 



■ i i | i i i | i i i | i i i | i i i | i i i | i i i | i i i | i i i | i i, 

0-40% centrality "; 




PHOBOS v 2 {EP} 
— 19.6 GeV 
-•- 62.4 GeV 
-•- 130 GeV 
-•- 200 GeV 



y 



-L. 



-L, 



i I ■ ■ ~l 



0.08 



-16 -14 -12 -10 -8 -6 -4 -2 

n-y 



1 — I I I I 1 1 1 1 1 — I I I I I I 




PHOBOS v,{EP> 
ALICE v 2 {2} 
ALICE v,{4> 



beam 



10 10 
\fs^ (GeV) 



Figure 15. a) Elliptic flow plotted versus r\ — j/beam for RHIC and the LHC, with 
for ybeam the y/s NN = 2.76 TeV value. The mesh shows the ry-range over which the 
elliptic flow is constant within statistical uncertainties, b) Elliptic flow versus beam 
energy. In both figures the uncertainties for ALICE are systematic uncertainties and 
the V2 in 0-40% centrality is obtained by averaging over events instead of over, the 
more commonly used, particle yield. 



extrapolation [59] from RHIC to the LHC was made (dashed line in Fig. 15 1) which 
predict an increase in v-i of ~50%, larger than predictions from most other models. 
PHOBOS measured v-i down to pt — using the eventplane method which at RHIC 
is similar to measuring t>2{2}. The measurements at rj — for the LHC are below 
the triangular extrapolation. However, the elliptic flow as function of pseudorapidity 
measured by ALICE in \rj\ < 0.8 (mesh in Fig. [TB^ i) is constant within uncertainties. If 
one takes into account that the 1*2(77) does saturate, like the multiplicity, at each energy 
around midrapidity then the longitudinal scaling might hold up to LHC energies. 

Figure [Hfo shows that the energy dependence of elliptic flow at midrapidity also 
seems to follow a rather simple scaling: the measured elliptic flow for four beam energies 
at RHIC shows a linear increase, extrapolating this to ^/s NN = 2.76 TeV results in a 
v 2 which agrees well with the ALICE measurement. 

The observed increase in the elliptic flow as function of beam energy is either due 
to an increase in ^-differential flow or due to an increase in the average transverse 
momentum of the charged particles. In most hydrodynamic model calculations the 
^-differential elliptic flow of charged particles does not change significantly j5U [52] . 
while the radial (azimuthally symmetric) flow does increase which leads to an increase 
in the average transverse momentum. The larger radial flow also leads to a decrease 
of the elliptic flow at low transverse momentum, which is most pronounced for heavy 
particles. Figure 16 1 compares the p%- differential elliptic flow of charged particles for 
three centralities at the LHC with STAR measurements at RHIC. It is seen that the 
^-differential elliptic flow is the same within experimental uncertainties. In Fig. [T6p 
the p t -differential elliptic flow measured at four beam energies is shown. The agreement 
of w 2(Pt) at these beam energies, which differ by almost two orders of magnitude, is 
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Figure 16. a) i>2{4}(pt) for various centralities compared to STAR measurements. 
The data points in the 20-30% centrality bin are shifted mp t for visibility, b) «2{4}(pt) 
for a collision energy range which covers almost two orders of magnitude |60j 



remarkable. Measurements of identified particle elliptic flow at these energies will reveal 
if this agreement can be understood in hydrodynamic model calculations. 

6. Summary 

In this review, I have shown that elliptic flow is one of the most informative observables 
in heavy-ion collisions. Nevertheless, the wealth of experimental information obtained 
from elliptic flow is far from being fully explored. The theoretical understanding of the 
experimental data is rapidly improving as is our understanding of the dynamics in heavy- 
ion collisions and the properties of the new state of matter, the quark gluon plasma. 
New high quality data from the LHC recently became available which shows that at 
LHC energies elliptic flow can be studied with unprecedented precision. This is because 
of the increase in particle multiplicity and also because of the increase in the flow signal 
itself. Due to the expected longer life time of the quark gluon plasma and the smaller 
contributions from the hadronic phase it is argued that the LHC is also better suited to 
determine rj/s of the partonic fluid (61]. Measurements of identified particle elliptic flow 
at the LHC and in particular the stronger mass dependence (splitting) of t>2(j>t) will 
be important to confirm the current theoretical picture. Additional constraints on 77/s 
can be obtained by measurements of the other anisotropic flow harmonics V3, V4 and v§. 
In the near future these measurements will become available and significantly increase 
our understanding of ultra-relativistic nuclear collisions and multi-particle production 
in general. 
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